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We report studies of controlled interactions of localized dissipative structures in a system described 
by the AC-driven damped nonlinear Schrodinger equation. Extensive numerical simulations reveal 
a diversity of interaction scenarios that are governed by the properties of the system driver. In 
our experiments, performed with a nonlinear optical Kerr resonator, the phase profile of the driver 
is used to induce interactions on demand. We observe both merging and annihilation of localized 
structures, i.e., interactions governed by the dissipative, out-of-equilibrium nature of the system. 


Localized structures coexisting with a homogeneous 
background are ubiquitous phenomena in extended dissi¬ 
pative systems driven far from equilibrium. These struc¬ 
tures consist of solitary excitations that manifest them¬ 
selves as electrical pulses in nerves [l|, concentration 
^ots in chemical reactions Q, oscillons in water waves 
[ 3 , 1^ and in granular matter Q, filaments in gas dis¬ 
charges [^, patches and fairy circles in vegetation H , or 
feedback and cavity solitons in nonlinear optics |8l-ll0l|. 
More generally, they are referred to as localized dissipa¬ 
tive structures (LDSs) or dissipative solitons [HQ. 


Like the conventional solitons of conservative inte- 
grable systems, LDSs can interact and collide with each 
other, sometimes with particle-like characteristics. But 
while conventional solitons always emerge unscathed 
from collisions 13|, LDSs can form bound states, merge 
into one, or even annihilate Q. These complex interac¬ 
tions arise from the non-integrability of nonlinear dissi¬ 
pative systems, and their study is of particular interest to 
better understand systems outside thermal equilibrium. 
Merging and annihilation of solitons have been exten¬ 
sively studied experimentally in non-integrable conser¬ 
vative systems, mostly with optical waves [Hill , but 
also, more recently, with matter waves [l6| . In contrast, 
although several authors have reported complex behav¬ 
iors of ensembles of LDSs in various settings, experi¬ 
mental observations have been uncontrolled and mostly 
qualitative (see, e.g., i El). It is only in gas dis¬ 
charges 171 and in vertically-driven fluids [J] that quan¬ 
titative measurements of the interaction laws have been 
obtained, with Q also resolving the merging dynamics. 
These two latter examples are realizations of, respec¬ 
tively, a reaction-diffusion system and a parametrically- 
driven damped nonlinear Schrodinger equation (NLSE) 
near the 2 :1 resonance. 


Here we report on a detailed numerical and experi¬ 
mental study of controlled merging and annihilation dy¬ 
namics of LDSs in a system described by an AG-driven 
damped NLSE near the 1:1 resonance. Experiments are 
performed in a nonlinear optical Kerr resonator, in which 
we can excite LDSs at selected and precise positions, and 


systematically induce their interactions. The interactions 
are triggered by manipulating the phase profile of the 
driver; the outcome controllably depends on the driving 
frequency and strength. Two LDSs either merge into 
one, or are both annihilated. In both cases, we tempo¬ 
rally resolve the collision dynamics and clearly observe 
the dissipative nature of the interaction through analysis 
of the energy balance. 

To better illustrate our experimental findings, we start 
our discussion by presenting numerical results. In dimen¬ 
sionless form, the AG-driven damped NLSE reads 

= + ( 1 ) 

This equation represents in our case the mean-field be¬ 
havior of a Kerr resonator [l^! but is also the small am¬ 
plitude limit of the AG-driven sine-Gordon equation Q. 
It has applications in non-equilibrium systems ran gin g 
from plasma physics to Josephson junctions jH |. 
highlighting the general applicability of our study. The 
equation can be cast into an autonomous form by sub¬ 
stituting = '0(a;, 

iipt + +'^xx =-itp + + iS (2) 

which will be used throughout this Letter. 

Depending on the driving strength S and its frequency 
A, Eq. (m exhibits a ran ge o f solutions, which have been 
extensively investigated jlSl . Q. Briefly, the simplest 
steady-state solutions are homogeneous {tpx = 0), and 
they satisfy the well-known cubic steady-state equation 
A = y3_2Ay+(A2-Hl)r with A = jS'p and F = |'0p. 
The steady-state curve {Y vs. A) is single-valued for 
A < v^, whereas for A > -x/S it assumes an S-shaped 
hysteresis cycle with three possible states. The latter 
range is of more relevance to our experimental configu¬ 
ration Q, and thus the focus of this Letter. Only two 
of the three states that exist for A > -x/S are homoge¬ 
neously stable (bistability): the negative slope branch is 
always unstable. For F > 1, the upper branch exhibits a 
Turing-pattern instability (also known as modulation in¬ 
stability) with respect to inhomogeneous perturbations, 
which can lead to the formation of a stationary periodic 
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Figure 1. (color online). Numerically simulated dynamics of 
induced LDS interactions for (j}{0) = 0.5 rad and xq = 30. In 
(a) [A, 5'o] = [2.91,1.87] and two LDSs merge into one; in (b) 
[A, So] = [3.64, 2.10] and two LDSs annihilate one another. 


pattern [l^ . LDSs can manifest themselves under condi¬ 
tions of coexistence of a patterned solution and a stable 
homogeneous solution (viz. Y < 1). They can be un¬ 
derstood to coincide with the patterned solution over a 
finite region in x, and with the homogeneous solution 
elsewhere (2^ . 


We are interested in the dynamics that take place when 
two LDSs collide. Unlike conservative solitons, widely- 
separated LDSs of the AC-driven NLSE are phase-locked 
to the driver, and thus all of them possess identical 
traits (for given X and A), including frequency (A) and 
velocity. Accordingly, unassisted collisions occur only 
when two LDSs are sufficiently close to interact attrac¬ 
tively 23|, yet such interactions are difficult to explore 
controllably. Inducing collisions by suitably modulat¬ 
ing the phase of the driver addresses that issue. 

Specifically, given S{x) = So exp[f(^(a;)], an LDS at xl 
will move towards the local maximum of (j}{x) with a 
drift velocity of dxn/dt = 4>'(xi,) A collision is 

thus observed when exciting, for example, two LDSs on 
opposite sides of a local maximum of 4>{x) [^ . 


To illustrate such induced collisions, we numerically 
integrate Eq. m using the split-step Fourier method. 
We assume a Gaussian driver phase prohle (j){x) = 
^(0) exp(— x^/xq). To create LDSs symmetrically dis¬ 
tributed about the phase maximum at x = 0, we use 
the initial condition -(/^(x, 0) = [sech[-\/A(x — xl)] + 
sech[\/A(x-|-XL)]] . The initial LDS separation 2 xl = 70 
is chosen to be much larger than their characteristic 


width (^ 1/A) B M [ 2 ^ so as to avoid any inter¬ 
actions during the transients leading to the LDS for¬ 
mation. Figures [T][a) and (b) show typical results for 
two different sets of driver frequency A and strength 
So, as listed in the caption, and with (/)(0) = 0.5 rad 
and xo = 30. These parameters are chosen to replicate 
our experiments. For clarity, the figures neglect the ini¬ 
tial portion of the simulation (which lasts for more than 
t = 1700), during which the two LDSs slowly approach 
each other from their initial separation of 2xl. In both 
cases, it can be seen that the LDSs drift towards each 
other until they are close enough to interact. The out¬ 
come of the collision is, however, markedly different. In¬ 
deed, for [A,S'o] = [2.91,1.87] the two LDS merge into 
one [Fig. [T][a)], whilst for [A, S'o] = [3.64,2.10] the intra¬ 
cavity field after the interaction is globally reduced to the 
homogeneous solution, i.e., the two LDSs annihilate one 
another [Fig. [T][b)]. 


It is apparent that the interactions depend on the pa¬ 
rameters of the driver. This has been numerically ex¬ 
plored further by systematically varying our four control 
parameters. A, So, (/(O), and Xo, over a wide range. We 
have found that A and So mainly govern the outcome 
of the collision, while the phase modulation parameters 
/>(0) and Xq mostly determine the speed at which the 
LDSs approach each other, i.e., set the timing of the col¬ 
lision. In Fig. [5J we summarize the observed outcome 
of the interaction as a function of So and A for the 
same driver phase modulation (/(x) as above. As can be 
seen, merging (green) and annihilation (blue) occurs in 
clearly distinct, but adjacent, regions. For a given driving 
strength So, the system favors annihilation over merging 
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Figure 2. (color online). Results from numerical simulations 
illustrating the outcome of LDS interactions as a function of 
driving frequency A and strength So for ((>(0) = 0.5 rad, and 
Xo = 30. Each solid dot represents a distinct simulation. No 
LDSs exist in the white-colorcoded area. 
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at higher driving frequencies. This can be related to the 
closer proximity to the folding point at Ac ^ tt^I 5 * 01 ^/ 8 . 
beyond which LDSs cease to exist in this system |19l| . 
Interestingly, in the area marked “bound states” no col¬ 
lision occurs. Instead, the two LDSs form a stable bound 
state repulsive interactions of the LDSs resist the 

drift induced by the driver phase modulation. Not sur¬ 
prisingly, this region slightly grows at the expense of the 
“merging” region when a shallower phase modulation is 
used (the merging/annihilation boundary is mostly unaf¬ 
fected). In the grey region, labeled “breathing”, the indi¬ 
vidual LDSs exhibit breathing as a result of an underlying 
Hopf bifurcation 0 . Their interaction can lead either 
to merging or annihilation, depending on the phase of 
their breathing at the onset of the collision. The bound- 
state and breather regimes will not be further discussed 
in this Letter because experimental limitations currently 
prevent us from observing them. 

We now describe our experimental configuration, im¬ 
plemented in the optical domain. Specifically, we induce 
controllable LDS interactions in a coherently-driven pas¬ 
sive optical fiber resonator that exhibits instantaneous 
Kerr nonlinearity. In the high-finesse limit, this system 
is known to be governed by Eq. 0 , with ■)/) representing 
the slowly-varying envelope of the electric field [l^ . The 
LDSs of such Kerr resonator have been observed exper¬ 
imentally before and are usually referred to as temporal 
cavity solitons [l^ . These are pulses of light that contin¬ 
uously circulate in the resonator, yet remain stationary 
in a reference frame that is moving at the group velocity 
of the driving light in the fiber. The transverse coordi¬ 
nate X in Eq. is thus a “fast-time” x ^ t that is de¬ 
fined in such a reference frame and that allows to describe 
the temporal profile of the field envelope. In contrast, t 
is a “slow-time” that describes changes in the field en¬ 
velope over consecutive roundtrips around the resonator. 
The normalization is such that dimensional time-scales t' 
and t' (units of s) and the electric field envelope E{t\ t') 
(units of W^/^) are related to the dimensionless variables 
in Eq. ([2]) by 0 


t' 

t = a —, 
tn 



i, = e^Ie, 

V a 


( 3 ) 


Here % is the roundtrip-time of the resonator, a is 
equal to half the percentage of total power loss per 
round-trip, L is the resonator length, and (< 0 ) and 
7 are, respectively, the anomalous group-velocity dis¬ 
persion and Kerr nonlinearity coefficients of the fiber. 
The driving strength Sq is related to the power Pin of 
the continuous-wave (cw) laser driving the resonator as 
S'o = (PinlLO/a^y^^, where 9 is the intensity transmis¬ 
sion coefficient of the coupler used to inject the field into 
the resonator. Finally, A characterizes the frequency de¬ 
tuning of the cw driving laser at ui from the closest res¬ 
onator resonance at wq, A ~ tR(wo — io)/a. 


Pattern generators 



Figure 3. (color online). Experimental setup. cw: 
continuous-wave, EDFA: erbium-doped fiber amplifier, BPF: 
band-pass-filter, SMF: single-mode fiber. 


A detailed schematic of our experimental setup is 
shown in Fig. [S] Overall, it is similar to the one used 
in [ 2 ^ . As a coherent driver, we use a narrow linewidth 
cw laser at 1550 nm wavelength, which is amplified up 
to 1.14 W using an erbium-doped fiber amplifier (EDFA) 
before being coupled into the resonator by a 90/10 fiber 
coupler {6 = 0.1). Noise accumulated during the ampli¬ 
fication stage is mostly removed with an optical band¬ 
pass filter (BPF). The resonator is composed of 100 m 
of standard silica single-mode fiber (SMF), with P 2 = 
—21.4 ps^/km and 7 = 1.2 W“^km“^. It also incor¬ 
porates an optical isolator to prevent resonance of stim¬ 
ulated Brillouin scattering radiation, and a 99/1 fiber 
coupler through which the intracavity dynamics are mon¬ 
itored with a fast photodiode and a real-time oscilloscope. 
The overall finesse of the resonator was measured to be 
F = 7 r/a ^ 21.5. The BPF at the 1 % output filters out 
the homogeneous cw background that coexists with the 
LDSs, thereby improving the signal-to-noise ratio of our 
data [l0|. The resonance frequencies of our optical fiber 
ring generally exhibit fluctuations due to environmental 
perturbations. To maintain a fixed A, we therefore ac¬ 
tively actuate the driving laser frequency to follow any 
changes in the resonances, by locking to a set level the 
optical power reflected off the resonator input 00. 
Changing the lock point allows us to controllably ad¬ 
just A, but we remark that the accuracy with which we 
can do so is insufficient to explore the formation of bound 
states since they manifest themselves over a narrow range 
of driver frequencies (see Fig. [2]). In this context, we also 
note that, with our current configuration, we are unable 
to reach power levels required to explore interactions of 







































I 


breathing LDSs. 

To controllably induce LDS interactions, we phase 
modulate the resonator driving field with a 10 GHz 
electro-optic modulator. The modulator is driven by 
one of two 10 GHz programmable pattern generators, se¬ 
lected with an electronic switch. The pattern generators 
are synchronized to each other by a single external clock, 
such that the repetition rate of their output patterns is 
identical to the resonator free-spectral range. The first 
generator (left in Fig. [3]) is configured to produce a pat¬ 
tern of two 130 ps full-width-at-half-maximum electronic 
pulses with 200 ps separation. These are fed to the phase 
modulator in the initial stage of the experiment. During 
that stage, we mechanically perturb the resonator, which 
results in the direct excitation of two LDSs at the two 
phase maxima [s^. After the LDSs are stably formed 
and trapped at the maxima we activate the elec¬ 
tronic switch and the phase modulator feed is abruptly 
changed (within a few nanoseconds) to the output of the 
second pattern generator. That generator is set to pro¬ 
duce a pattern made up of a single pulse whose delay 
is adjusted to lie halfway between the two pulses gener¬ 
ated by the first generator. Accordingly, the LDSs in the 
resonator find themselves approximately symmetrically 
positioned about the new single maximum of the phase 
profile. As in the simulations of Fig. [U the LDSs thus 
start drifting towards that maximum, interacting once 
sufficiently close to each other. Note that the new driver 
phase profile takes a few photon lifetimes (^ 1 fis) to get 
imprinted inside the resonator after the switch, but that 
transient is negligible in regards of the interaction time. 

The temporal intensity profile of the intracavity light 
measured at the 1 % output of the resonator is recorded 
by the oscilloscope (triggered by the pattern genera¬ 
tors) as we operate the switch. Successive recordings 
are vertically concatenated and shown as density plots 
in Figs. Sla) and (b), which have been obtained for 
[A, S'o] = [2.91,1.87] and [3.64,2.10], respectively. The 
first 10 s of the measurements are very similar: two LDSs 
with 200 ps separation are stably trapped at the maxima 
of the phase pulses defined by the first pattern generator. 
After switching to the single phase pulse pattern (which 
occurs at t ~ 10 s) and inducing the LDS interaction, 
a single LDS is seen to remain for [A, S'o] = [2.91,1.87] 
whilst both disappear when [A,So] = [3.64,2.10]. These 
results are indicative of merging and annihilation, re¬ 
spectively, which is in agreement with numerical sim¬ 
ulations. Indeed, the simulation results in Fig. [T] use 
the very same parameters as the experiments here. Yet, 
these results are limited by the slow 1 frame/s acquisi¬ 
tion rate of the oscilloscope which does not reveal the 
transient energy balance dynamics. We have thus also 
recorded the roundtrip-to-roundtrip evolution of the in¬ 
tracavity energy on the real-time oscilloscope. Experi¬ 
mental results for merging and annihilation are shown 
as red circles in Figs.SKc) and (d), respectively, superim- 
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Figure 4. (color online), (a, b) Experimental density plots 
showing the evolution of the intracavity temporal intensity 
profile as two LDSs (a) merge into one, and (b) annihilate each 
other. The successive traces are recorded with a real-time 
oscilloscope at 1 frame/s. (c, d) The roundtrip-to-roundtrip 
evolution of the total intracavity energy during LDS (c) merg¬ 
ing and (d) annihilation. The experimental data is repre¬ 
sented by red circles, and results from numerical simulations 
are shown as blue solid lines. 


posed with results from numerical simulations (blue solid 
lines). Here we only show the portion of the dynamics 
corresponding to the final stages of the interaction (the 
energy stays approximately constant during the slow ap¬ 
proach of the two LDSs) and we have normalized the 
energy such that a single isolated LDS carries the en¬ 
ergy 1 a.u. We have also post-processed the numerical 
simulation results, extracted from data shown in Fig. [TJ 
to take into account the BPF at the resonator output 
as well as the limited bandwidth of our photodetector. 
The results in Fig. UJc, d) clearly confirm that merging 
and annihilation occur in our experiment and more im¬ 
portantly the dissipative nature of the interactions. For 
[A, S'o] = [2.91,1.87] the energy falls from two to one, im¬ 
plying merging; for [A,So] = [3.64,2.10] the energy falls 
from two to zero, implying annihilation. 

These results represent, to the best of our knowledge, 
the first realization of controllable interactions of local¬ 
ized dissipative structures. Our study also provides the 
first quantitative analysis of such interactions in an AC- 
driven nonlinear Schrbdinger system, and more generally, 
in any nonlinear dissipative system near the 1:1 reso¬ 
nance. We have numerically identified a diversity of in¬ 
teraction scenarios for different parameters of the system 
driver. Experiments performed in an optical resonator 































5 


show unequivocal evidence of possible selection of LDS 
interaction by the operator from merging to annihilation. 
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foundation. 
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